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Abstract 

In an effort to generalize the self-consistent Ornstein-Zernike approximation (SCOZA) 
- an accurate liquid-state theory that has been restricted so far to hard-core systems- 
to arbitrary soft-core systems we study a combination of SCOZA with a recently devel- 
oped perturbation theory. The latter was constructed by Ben-Amotz and Stell [J. Phys. 
Chem. B 108,6877-6882 (2004)] as a reformulation of the Week-Chandler- Andersen per- 
turbation theory directly in terms of an arbitrary hard-sphere reference system. We 
investigate the accuracy of the combined approach for the Lennard-Jones fluid by com- 
parison with simulation data and pure perturbation theory predictions and determine 
the dependence of the thermodynamic properties and the phase behavior on the choice 
of the effective hard-core diameter of the reference system. 



I. INTRODUCTION 



In its most general 



mation (SCOZA, can be obtained from a chosen liquid state theory by the 

introduction of an adjustable parameter (such as an effective temperature) and the 
subsequent imposition of consistency between two different routes to thermody- 
namics, typically the energy and compressibility routes. Simple as this prescription 
is, it has nevertheless proved to be highly effective and to lead to very accurate 
predictions for structure and thermodynamics throughout the temperature-density 
plane, and for the critical point and phase coexistence in particular Q, 0, 0, Q|. 
Indeed, SCOZA is one of only very few liquid state theories that do not develop 
serious problems in the critical region and even exhibit some form of scaling with 
non-classical, partly Ising-like critical exponents jjj. 

While this is by no means inherent to the concept of SCOZA, applications to 
continuum fluids have been restricted in a number of ways for both historical and 
practical reasons: Originally the theory grew out of the semi-analytic solution 
of the Mean Spherical Approximation (MSA) for hard core Yukawa potentials. 
Correspondingly, SCOZA has so far only been used with an MSA-like closure and 
with potentials composed of a hard core of diameter <th and an attractive tail 
va{t), a situation we will refer to as the "hard- attractive" (HA)one. As far as 
v a is concerned, the original restriction to a single Yukawa term [6] has gradually 
been relaxed, first by expanding the class of admissible tails to superpositions of 
two [^J or more Yukawa terms |lf| . 



ormulation, the Self- Consistent Ornstein-Zernike Approxi- 
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and then to those of Sogami-Ise form 



EL 



zv (r— <th) 



V 



for which semi-analytic results are also available, only to be fully overcome through 



tion 12]. Approximation of v&(r) as the superposition of a rather small number 
of Yukawa or Sogami-Ise tails — a non-trivial step for some of the artificial model 
potentials used in liquid state theory such as, e. g., the square well one — is thus 
no longer necessary, and it becomes possible to use exactly the same interaction 
as the simulations one may want to compare with. Not surprisingly, the flexibility 
so gained comes at substantial computational cost in solving the SCOZA partial 
differential equation (PDE). 

Despite this liberation of the form of the attractive tail, application of SCOZA 
has so far remained tied to a hard core reference fluid (marked by subscript H) 
with potential 



Only in the case of a bounded interaction like, e. g., the Gaussian one, v(r) oc 
exp(-a^), has the reference system effectively been eliminated by letting 
go to zero formally [13]. But this is not an option for simple liquids where the 
potential diverges for r — > 0, such as the classic Lennard- Jones (LJ) one, 



For this or more specific systems the infinitely strong hard core repulsion vn(r) 
present in SCOZA is hardly realistic. It is therefore desirable to generalize the 
theory to largely arbitrary, strong but soft short-range repulsion vs(r) to which a 



the introduction of a fully numerical solution of the Ornstein-Zernike (OZ) equa- 




(1) 




(2) 
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weak attractive tail v^(r) is again added to obtain the full "soft-attractive" (SA) 
system, 

^sa(>) = v(r) = v s (r) + v A (r). (3) 
Just as in an approximate SA-SCOZA recently proposed that relies on the 



virial theorem to gauge t 



le amplitude of the pair distribution function (PDF) at 



the molecular surface ^|, in the present contribution we also opt for the time- 
honored strategy of representing the soft repulsive core by hard spheres of state- 
dependent diameter. We therefore need explicit prescriptions for both the state- 
dependence of this diameter and for any correction terms that may be needed to 
account for the softness of the reference fluid. In doing so there is considerable 
atitude, and many different schemes have been proposed and used in the past 
I5I llfij ]. For example, the highly successful first order perturbation theory due 
to Weeks. Chandler and Anderson (WCA. 17]) solves both of these problems by 
determining <tr as a function of the temperature T = ks/P and the particle density 
p such that the Helmholtz free energies of the soft core reference system and the 
hard spheres coincide. Determination of the diameter is then rather involved and 
leads to a dependence not only on temperature but also on density. Even though it 
has been argued that p dependent cxr is fundamentally more appropriate 
practical advantages of a simple prescription that depends on temperature only 
often outweighs a slight loss in accuracy. Presumably this accounts for much of 



1UUJ 

18], 



the 



the lasting popularity of the Barker-Henderson (BH) theory according to which 
or(/3) is computed as j3] 



<7 BH (/3) = / (1 - e-^M) dr . (4) 
Jo 

In this context a recent first order perturbation theory due to Ben-Amotz and 
Stell (BAS, |3 Q) i s °f particular interest: While WCA consider a general ref- 
erence system that is in turn mapped onto an effective hard sphere fluid to which 
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attractions are then added, BAS start from a pure hard sphere reference fluid in- 
stead and add two terms to the free energy to account for v$ and v a, respectively. 
This hard core WCA, or BAS, theory has been shown to yield results of an accu- 
racy comparable to that of the original WCA theory while being insensitive to the 
precise choice of <th- This makes it possible to use an even simpler temperature 
dependence of the effective diameter o"h than that of Eq. (jlj) without loss of accu- 
racy [lij , and the resulting scheme combines the advantages of the WCA and BH 
theories. 

It is this BAS approach to the treatment of the reference system that we here 
set out to combine with SCOZA along the lines of a suggestion first put forward 
by Raineri and co-workers j^J. After a short reminder of SCOZA and BAS theory 
and a presentation of their combination (section |H} we apply our BAS SCOZA 
to the LJ interaction (j2J), comparing our results with simulation data as well as 
the predictions of BAS theory ( section llll|) . It turns out that we can well repro- 
duce earlier BAS results and, in fact, improve upon them in slightly supercritical 
isotherms. As far as phase coexistence is concerned, however, the residual diam- 
eter dependence of the free energy translates into an uncertainty of the pressure 
P and chemical potential [i that renders the coexistence curve more sensitive to 
variations in cth than expected. As shown in the appendix, this is related to the 
divergence of the compressibility at the critical point. Without a criterion to fix 
<Jk{T) unambiguously, and in the face of further problems connected to the spin- 
odal of the underlying SCOZA computation if the hard-sphere diameter of the 
reference system gets too small, the BAS SCOZA considered here is thus found to 
lack the faculty of predicting the L J phase diagram with the accuracy one has come 
to expect from SCOZA. In our concluding remarks we consider the prospects for a 
SCOZA adapted to non-hard sphere reference systems more generally ( section llVjl . 
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II. BAS CORRECTIONS TO SCOZA 



The starting point of BAS theory, and hence also of our combined BAS SCOZA, 
is the splitting (J3J) of the full potential into its repulsive and attractive components. 
In the present work we adopt the conventional WCA prescription of separating 
them at the position r min of the potential's minimum, i. e., 



v s [r 



v A (r) 



v(r) - v(r min ) for r < r min 
otherwise 

(5) 

) for r < r min 
v(r) otherwise; 



in the LJ case, r min = \/2a and v (r min ) = — e. The Helmholtz free energy of the 
fluid is customarily written as the sum of the ideal gas contribution A ld and the 
excess free energy A ex , 

A = A id + A cx . 

In computing v4 ex , BAS start from a hard sphere reference system and split v4 cx 
into the hard sphere term Aff and two correction terms corresponding to Vs and 
va, respectively, 

A ex = A^ + A A s + AA A . 

The term AA§ implementing the difference between hard (H) and soft (S) 
repulsion is given by [lq 

(3AA S o [°° . . . , . . a 



N 



2vrp / (gu(r) - gs(r)) r dr, (6) 



where the g(r) denote the respective PDFs. As in WCA and BAS theories, we re- 
late g$ to <7h by the assumption of equal cavity function y(r) = g(r) exp (+{3v(r)). 
This is expected to be quite harmless for the rather stiff LJ case considered here 



6 



2l\ . We have checked that the three-term approximation for yu{ r ) inside the 
core used in Ref. [15l gives essentially the same results as a four-term one due 

n 

to Grundke and Henderson \22\ . Outside the core we adopt Waisman's prescrip- 
tion for the hard sphere structure functions: The direct correlation function ca(r), 
r > oh, is approximated by a single Yukawa tail with parameters chosen so as to re- 
produce the Carnahan Starling equation of state through both the compressibility 
and the virial routes j3] , cf. appendix A of Ref. (f| . 

The remaining contributions to the free energy are and For the 

former one may use the well-known Carnahan-Starling resul t | 24l| . whereas the 
latter is given in BAS theory as an integral over gu(r) va(t) (la|. On the other 
hand, A^ A = A^ + AAa corresponds to the superposition of hard core repulsion 
and long-range attraction, i. e., to the potential 

VnA{r) = v n (r) +v A (r). 

This is precisely the HA situation SCOZA has long been applied to with excellent 
results. In combining BAS theory and SCOZA it is therefore natural to use SCOZA 
for the description of this HA sub-problem, and for the computation of A# A in 
particular. By adding AA$ we then arrive at the full excess free energy, from 
which any remaining thermodynamic quantities may be derived. This recipe was 
first proposed in Ref. j^. 

The first step in computing the thermodynamics of the full SA system along 
an isotherm at inverse temperature (3 within the combined BAS SCOZA is the 
determination of the effective hard sphere diameter cth- In view of the purported 
insensitivity of BAS theory results to the precise choice of <th, we generally employ 
a simple Boltzmann factor criterion (BFC) to fix the temperature dependence of 
<7h [3, According to this criterion, we choose oh(/3) so as to keep /3i>s(o"h) 
constant, exp(— /3ws(c"h)) = 1/obfc- Here we have introduced a parameter obfc 
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whose variation allows us to change <th in a consistent way for all temperatures. 
In Ref. ^(|, BAS found virtually unchanged results for the parameter range 2 < 
obfc < 5, which includes the a priori preferred value cibfc = e - In some 

of the computations we also use the BH prescription (J3J. 

With <th so determined, the SCOZA part can be handled in the usual way Q| 
except that the WCA splitting (jSj) effectively precludes a multi- Yukawa or Sogami- 
Ise representation of v A (r) for r > o"h- We thus have no choice but to turn to an 
implementation where the OZ relation, along with the customary SCOZA closure 



is solved numerically [l2j. The first of the above relations expresses the impenetra- 
bility of the effective hard cores whereas the second one describes the effect of the 
attractive tail on the HA direct correlation function cha(^) outside the core beyond 
that of hard spheres, cn(r). Comparison immediately shows that this corresponds 
to the MSA closure with j3 replaced by a state dependent effective inverse tem- 
perature K((3,p). The latter is fixed through the requirement of thermodynamic 
self-consistency as embodied in the SCOZA PDE 



Here 1/xua = 1 — P La cnA( r )d 3 r so that k = /?Xha/p is the isothermal com- 
pressibility of the HA system as evaluated from the compressibility route, and 
m H a = U^ A /V = 2np 2 gnA{r)vuA{r)r 2 dr is the energy route result for the ex- 
cess, or configurational, internal energy per unit volume. Solution of this PDE 
subject to suitable boundary conditions at infinite temperature as well as at van- 
ishing and high density yields both the structural and the thermodynamic proper- 
ties of the HA system throughout the domain of the PDE. Only below the critical 



guA(r) =0 for r < (j H 

chaO) = c H (r) -Kv A (r) for r > a H , 



(7) 
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temperature is there a region of instability or non- convergence that must be ex- 
cluded from the integration. It is customary to do so through the imposition of 
an additional boundary at the spinodal so that no SCOZA results are available 
inside the HA spinodal. For a detai 
refer the reader to the literature 



ed description of the numerical procedure we 



mm 

For use in BAS SCOZA, integration of the SCOZA PDE is stopped upon reach- 
ing the temperature corresponding to <th, and internal energy, free energy, pressure, 
and chemical potential of the HA system are extracted along this isotherm. Fi- 
nally, we add the BAS correction terms following from AA$ in order to arrive at 
the final SA results for these quantities. 

Determination of the SA phase diagram involves repeating this procedure for all 
the diameters corresponding to the temperature range of interest. As both the HA 
free energy so obtained and the BAS correction term AA$ are unique functions of 
temperature and density, so is their sum A ex , the excess free energy of the full SA 
system. The SCOZA PDE is then trivially fullfilled as long as the thermodynamic 
quantities are computed by differentiation of A ex . On the other hand, the structure 
of the SA system differs from that of the HA system computed with SCOZA and 
is therefore not accessible. Consequently, neither the compressibility route nor 
the energy route to thermodynamics can be evaluated for the SA case within the 
present approach. In this sense, BAS SCOZA must not be seen as a soft core 
version of SCOZA but rather as a BAS theory built upon a SCOZA foundation. 

III. APPLICATION TO THE LJ FLUID 

In order to gauge the performance of BAS SCOZA as described in the preceding 
section we now turn to its application to the LJ potential (J2J) with the WCA 



splitting (j3J). For comparison purposes we make use o 
simulation data published by Johnson and co-workers 



the molecular dynamics 
and by Lotfi and co- 
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workers [26] as well as Monte Carlo (MC) results due to Wilding [27]. While 
the former two directly relate to the full interaction v^j of Eq. (|2*|). the latter are 
for the cut but not shifted potential t>Lj(V)©(2.5c7 — r) without corrections (0 is 
Heaviside's function). We have therefore exploited the flexibility brought about 
by the fully numerical solution of the OZ relation and performed BAS SCOZA 
computations with both of these potentials, depending on the simulation data we 
want to compare with. 

In presenting their theory 



as well as in their comparative study of various 



thermodynamic perturbation theories BAS focussed on high densities, and 
mostly on rather high temperatures: With only a single subcritical temperature 
and with densities that are at least twice the critical one, all of the four states 
repeats cohered in Re, Q ar e q nite far from the critical point and the 
spinodal. Presumably they were chosen because it is in dense or hot systems 
that particles explore distances around <th effectively so that high temperatures 
and densities present the most challenging test for a perturbative description of 
soft cores. The effects of variations of <th generally increase with density for the 
isotherms studied in Ref. Q] as well as in our own calculations. 

For a comparison of BAS SCOZA with pure BAS theory (or actually, AA$ be- 
ing identical in both cases, of their HA parts only) these states are only of limited 
interest: At high density, the pair structure is dominated by packing effects so that 
any difference between (3 and the parameter K of the closure ((Zj) hardly affects the 
HA energy integral and the free energy obtained from it by thermodynamic inte- 
gration. For the lower density states, on the other hand, the high temperatures and 
great separations from the critical region imply that the SCOZA self-consistency 
problem can hardly have changed Kj (3 from unity appreciably and so again render 
the HA structure and thermodynamics of BAS SCOZA equivalent to pure BAS 
theory. Any remaining discrepancies in the thermodynamics must be attributed 
to the difference between first-order perturbation theory and what is essentially 
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the MSA for the HA problem. It is therefore no surprise that BAS SCOZA repro- 
duces the BAS results for these states very well. A detailed comparison of the data 
underlying Fig. with the corresponding Fig. 2 of Ref. does, however, seem 
to hint at a sensitivity of A ex to variation of oh that is slightly reduced vis-a-vis 
BAS theory at T* = ksT/e = 2.81, p* = per 3 = 0.85, and slightly increased for the 
smallest <th values considered at T* = 0.75, p* = 0.84; for the latter temperature, 
however, oh ~ 0.94 corresponds to obfc ~ 10 2 , which is far outside the normal 
range, 2 < a BFC < 5. For the two remaining states, T* = 3.05, p* = 1.1 and 
T* = 1.35, p* = 0.65, any differences are too small to be made out from Fig. 2 of 
Ref. Q- 

SCOZA differs from MSA mainly in the critical region and in the vicinity of the 
spinodal where K strongly deviates from /3. In order to see a genuine SCOZA con- 
tribution, as opposed to merely gauging the accuracy of a first order perturbation 
theory for the HA subproblem, our interest is naturally drawn closer to the HA 
critical point, and thus also closer to the critical point of the full SA system. Of 
the isotherms displayed in Fig. 4 of Ref. [lfj], the one at T* = 1.35 should already 
be sufficiently close to the critical temperature of the full LJ interaction that has 
been estimated as T* = 1.310 (26 1 a nd 1.313 by molecular dynamics, and 
T* = 1.3120(7) [2J and 1.326(2) [2g| by Monte Carlo methods. Unfortunately, 
BAS only show Z = f3P/ p in the figure so that it is difficult to discern whether 
a van der Waals loop is present, nor do they address this question directly. On 
the other hand, the standard by which the performance of BAS theory is judged 
in Ref. ^| is the classic WCA theory. As all of BH theory, WCA theory, and a 
thermodynamically self-consistent variant of the latter due to Lado display almost 
identical van der Waals loops at even higher temperature (T* = 1.4, cf. inset in 
Fig. 5 of Ref. [3|), we may safely assume that this is true for the BAS result, too. 
As can be seen from the inset in our corresponding Fig. the same isotherm is 
correctly predicted to lie above the critical temperature in BAS SCOZA and is in 
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essentially perfect aggreement with the simulation data for 2.3 < obfc — e - 

By way of contrast, the same diameters give substantial deviations and show a 
marked trend towards a van der Waals loop for smaller obfc when BAS SCOZA 
is replaced by "BAS MSA", i. e., if K is restricted to coincide with [3, cf. inset 
in Fig. El While an adjustment of obfc so as to achieve good agreement with 
the simulation data in this part of the phase diagram is certainly possible, the 
resulting diameter may be too small to describe the packing effects at high density 
correctly. At any rate, for slightly supercritical isotherms like that at T* = 1.4, 
BAS SCOZA not only leads to a better agreement with the pressure data but also 
proves less sensitive to variations of the diameter. It is therefore superior not only 
to pure BAS theory (which might be explained by the qualitative difference of per- 
turbation theory vs. integral equations) but also to MSA with BAS corrections. 
The improvements must therefore be attributed to SCOZA's self-consistency re- 
quirement. 

Returning to BAS SCOZA, even for this narrow obfc range there appears a 
systematic, if small, pressure difference that increases with density. Although it is 
not obvious how to relate the sensitivities to variation of <7h of different quantities, 
comparison of Figs.|2]and|n]seems to indicate that the pressure depends on o"h more 
strongly than the free energy. A similar conclusion can be drawn for the chemical 
potential by comparing Figs. ^ and EJ This greater sensitivity of P and fi relative 
to A ex for most states, and for comparatively low temperatures in particular, is 
also demonstrated in Tab. |J where <th induced changes are related both to the 
absolute value and to the density dependence of the respective quantities. It can 
also be seen in pure BAS theory: For the lowest isotherm in Fig. 4 of Ref. [l(| . 
T* = 0.74, a slight variation of <Th from cr^p^ by only two per cent is sufficient 
even to change Z from a convex to a concave function of density for p* w 1. 

Proceeding to the critical region and phase coexistence, in BAS SCOZA the 
binodal is located as in pure SCOZA viz., by a search for densities of equal pressure 
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P and chemical potential /i; the critical point is identified with the locus where 
the gas and liquid branches of the binodal meet. As expected for a variant of 
BAS theory even if it makes use of SCOZA input, the critical behavior is not 
compatible with the Ising universality class: In the accessible temperature range 
the coexistence curve is not described well by the usual scaling form, and for the 
highest sub-critical temperatures the effective exponent (3 tends to values far larger 
than SCOZA's 7/20 [8], which in turn slightly exceeds the correct Ising value. This 
is also conspicuous from the forms of the binodals obtained from the BFC diameter 
cbfc with several values of czbfc as well as the BH diameter <jbh as displayed in 
Figs. El and El for the truncated and the full L J potentials, respectively. 

Figs. El and El also bring out the gravity of the diameter sensitivity of P and 
\i for the description of phase coexistence: There is a substantial variation in the 
locations of the upper parts of the coexistence curves when o"h is changed from 
a BF~c 2 ' (Fig- El) or o^yc^ 1 (Fig- EJ- It should be noted that these optimal values 
of obfc cannot be computed from the theory itself but are merely the results of 
comparisons with the simulation data for the coexisting densities for the truncated 
and full potentials, respectively. 

The particularly great sensitivity of the binodal at the highest temperatures 
can in fact easily be understood by linking the shift in the coexisting densities 
induced by a change in the effective hard core diameter o"h to the isothermal com- 
pressibility k at phase coexistence. Clearly, k diverges at the critical point and 
then decreases along the two branches of the binodal as we proceed to lower tem- 
peratures. Informally speaking, a higher compressibility means that the physical 
system must be compressed or expanded, and therefore its density changed, by a 
larger amount in order to offset a small change in pressure and chemical potential, 
which explains the more pronounced effect at higher temperatures. At the same 
time, the thermodynamic quantities are generally more sensitive to a change in 
<7h at higher densities, which explains the asymmetry of the effect between the 
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low and high density branches of the binodal. A more formal exposition of this 
reasoning can be found in the appendix. 

A different perspective on the strong <th dependence of the binodal is offered by 
the realization that an increase in <th renders the HA subproblem more strongly 
repulsive, thus lowering the HA critical temperature. In an exact theory, the cor- 
rection term AA$ strictly compensates this shift of the critical point. From a 
first-order perturbation theory, however, we can expect only a partial compensa- 
tion, as is clearly demonstrated in Fig. [7|where we separate the SCOZA and BAS 
contributions to the free energy. Consequently, a change in obfc (or, more gener- 
ally, in the prescription for crn(P)) still leads to a systematic shift of the SA critical 
temperature. For any given isotherm, this necessarily corresponds to changes of 
the coexisting densities that are most pronounced close to the critical point where 
the binodal gains in width most rapidly. 

Such shifts of (3 C can actually be inferred from Figs. |S] and HI The reason they 
cannot be seen directly is connected to the specific way in which BAS SCOZA 
combines the theories it is built upon, viz., by adding a BAS correction term to 
a SCOZA free energy: As mentioned in the introduction, SCOZA suffers from 
a region where no solution can be found due to either stability or convergence 
problems. This region is buried within the HA spinodal and is customarily elim- 
inated through the imposition of an artificial spinodal boundary condition. No 
BAS SCOZA results can therefore be obtained inside the HA spinodal. Depend- 
ing on the diameter, this hole in the solution may lie well within the SA binodal. 
For czbfc above a certain threshold, however, or more generally whenever <th(/3 c ) is 
too small, the HA interaction is not sufficiently repulsive so that the HA spinodal 
rises above the SA binodal, the upper part of which is then no longer accessible. 
As can be seen from Figs. El and this is the case for those diameters that come 
close to the simulated phase diagrams. The hole in the solution for temperatures 
below the HA critical one can also be seen from the missing intervals in the lowest 
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temperature isotherms displayed in Figs. El El and [TU1 

Let us now consider the configurational internal energy U' 



cx 



(d(3A™/d(3) p . 



With A ex = A^ A + AAg we obtain U ex as the sum of the HA internal energy 
£/£\ extracted directly from the SCOZA computation with fixed diameter and the 
correction AUs implementing the difference between hard and soft repulsion. The 
latter is composed of the temperature derivative of the BAS correction AA$ and a 
term related to the temperature dependence of Af[ A through <th (and thus to the 
steepness of v$ at the effective molecular surface), 



Here, the derivative dun/d/? is to be evaluated according to the convention used 
to fix <Th(/3), %■ e., from Eq. (jU) for <TbH) or at constant a BFC for cxbfc- 

From Fig. IHlwe see that the overall diameter dependence of the excess internal 
energy is qualitatively similar to that of the excess free energy, cf. Fig. ^ Not 
surprisingly, it is most pronounced at the highest densities whereas the temper- 
ature is of minor importance for the states considered in Fig. |HJ This can easily 
be understood in terms of packing effects that dominate the pair structure. In 
Fig. El we again separate the SCOZA and BAS contributions to the internal en- 
ergy. Interestingly, the latter compensates the diameter dependence of the former 
not even approximately as was the case for the free energy (Fig. HJ). Instead, AUs 
is essentially constant over the whole cth range shown, and the variation of U ex is 
due to Uft A alone. As Fig. El shows, the particularly low sensitivity of U cx with 
respect to cth m the upper panel of Fig. El is a mere consequence of the transition 
between density ranges where U ex rises or falls, respectively, with growing <7h- In 
the lower panel, the internal energy again varies by several per cent in the <7h range 
displayed. 



ex 
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Both pressure P and internal energy f/ ex as computed within BAS SCOZA with 
2.3 < obfc < e are generally in good agreement with the simulation results for the 
sample isotherms displayed in Figs. |2] and EH Only a detailed comparison reveals 
some systematic deviations that seem to hint at a moderate state dependence of 
the optimal value of obfc : At high densities and supercritical temperatures U ex is 
best reproduced with czbfc — e ; an d the simulated pressures indicate an obfc even 
slightly larger than e. By way of contrast, phase diagram (Fig.EJ), pressure (Fig.|2J) 
and internal energy (Fig. ITUj) all agree that <2bfc = 2.3 is the optimum value for 
the full LJ interaction for states for which SCOZA's self-consistency requirement 
is expected to be relevant, i. e., for T < T c as well as for T ~ T c , p ~ p c . A similar 
level of agreement is also expected for the truncated potential with a slightly larger 
diameter, Qbfc — 2.2, cf. Fig. |3J This difference of about 0.1 per cent in o"h seems 



remarkably smal 
a cut may have 



_given the pronounced influence on the critical parameters such 



301 ] . and in particular given the respective critical temperatures 



T* < 1.2 for the cut potential [22] as opposed to T* > 1.3 for the full interaction 
|2a bd, ba I23 , c/. Figs. E and El 



IV. CONCLUSION AND PERSPECTIVES 

As we argued at the end of section HH BAS SCOZA should be regarded as a 
SCOZA-based variant of BAS theory rather than as a modified SCOZA. From this 
point of view our expectations regarding the performance of this combined theory 
are well fulfilled: Good agreement with simulation results for binodal, pressure 
and internal energy can be achieved with a judicious choice of the parameter gibfc; 
we find a significant improvement over pure BAS theory at slightly supercritical 
isotherms; and away from the critical region and phase coexistence both schemes 
give essentially equivalent results. Unfortunately, the transition from pure BAS 
theory to BAS SCOZA entails a dramatic increase in computational cost as the 
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solution of a non-linear PDE in (3 and p is required just for the results along a 
single isotherm. For the computation of a full phase diagram, BAS SCOZA theory 
is far more demanding than even fixed-diameter SCOZA where a single integration 
of the PDE yields the results at all densities and temperatures. 

On the other hand, the present scheme was originally proposed as a soft-core 
extension of SCOZA rather than as an improved BAS theory [^(j]. And indeed, 
BAS SCOZA is able to handle SA fluids using the true interaction whereas previous 
applications of conventional SCOZA to the LJ fluid had to rely on heuristic, and 
by no means small, corrections to form and amplitude of the potential outside the 
fixed hard core [^J yielding a critical temperature that is too low by about five per 
cent Still, the huge increase in computational cost, the problems connected to 
the HA spinodal, and the loss of all structural information may seem an inordinate 
prize to pay for this advantage of BAS SCOZA. 

Most importantly, though, the distinctive feature of SCOZA, viz., the thermo- 
dynamic consistency requirement implemented by an effective temperature that 
replaces the true temperature in the MSA closure for the HA system, shows up in 
the results only in the proximity of the critical point. It is there that we have to 
look for SCOZA-specific improvements: In the rest of the (/3, p) plane SCOZA is 
essentially equivalent to MSA, and any improvements over pure BAS theory merely 
reflect the superiority of integral equations over perturbation theory in evaluating 
the difference AAa of the free energies of the H and HA systems. Close to the 
critical point, however, the chief advantage of BAS theory, viz., the near- constancy 
of the thermodynamic results under variation of cth breaks down, and the BAS 
SCOZA binodal strongly depends on o"h(/3). The value of the parameter obfc 
thus greatly influences the predicted phase behavior; at the same time, it cannot 
)e computed from the theory itself, and the a priori attractive choice of obfc = e 
2fl| is certainly unsatisfactory in the critical region. Further considering that the 
BAS SCOZA binodals are far from compatible with the Ising universality class 
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we therefore conclude that the theory as presented here cannot be used to predict 
the critical properties and liquid-vapor phase behavior from first principles with 
anywhere near the accuracy usually associated with SCOZA. 

As for the feasability of a BAS-based SCOZA rather than a SCOZA-based 
BAS theory, the main obstacle to using BAS results to "provide reference-system 
input" for SCOZA is the purely thermodynamic nature of the BAS term A As 
that does not allow us to relate the structures of the HA and SA fluids. Without 
.fT-dependent SA structures, however, neither the energy nor the compressibility 
routes to the SA thermodynamics can be evaluated, and the question of their 
consistency becomes meaningless. Indeed, inclusion of the soft core contribution 
in the thermodynamic self-consistency problem is likely to be of prime importance 
for the performance of an eventual soft core variant of SCOZA 14 1. 

What then, one may ask, are more promising ways towards a SCOZA capa- 
ble of describing SA systems? For this we see two options: One possibility is to 
consider formulations where the introduction of an effective hard core reference 
system still allows the analytical short cuts for the solution of the OZ equation 
available for the multi- Yukawa or Sogami-Ise fluids to be used and that include the 
softness in the self-consistency problem, if only approximately; one such approach 
has recently been proposed [14J. On the other hand, if the hard- won freedom to 
choose almost arbitrary functions v a and i>g for the attractive and soft-repulsive 
parts of the interaction is to be retained, the OZ relation must be solved numer- 
ically at any rate. In this case there is no need for a hard core reference system, 
and non-perturbative ways of treating the soft repulsive reference fluid may be 
more appropriate and less wasteful in terms of computing resources and structural 
information. We intend to pursue this topic in a further study. 



18 



Acknowledgments 

AR wishes to thank Johan H0ye for many stimulating discussions and grate- 
fully acknowledges financial support from the Austrian Science Fund (FWF) under 
project J2380-N08. This work was also supported by the Austrian Science Fund 
(FWF) under Project No. P17178-N02. 

APPENDIX A: DIAMETER SENSITIVITY OF THE BINODAL 

Let us consider a subcritical isotherm, (5 > f3 c , with vapor and liquid coexisting 
at densities p v and p h respectively. Along this isotherm the exact free energy is 
a function of density only, A = A(p). The BAS free energy, on the other hand, 
is computed in an approximate way and makes use of an effective hard sphere 
reference system with density independent diameter or- It therefore acquires an 
additional dependence. To obtain the free energy at slightly different values of ctr 
and p we can use the first order Taylor expansion 

A(p + Ap, a H + Acr H ) = A(p) + Aa H A(p) + Ap A\p) + . . . . 

In this appendix we use dots to denote differentiation with respect to oh and 
primes for differentiation with respect to p. For the benefit of a compact notation 
we will also use subscripts v and / to indicate functions evaluated at p v and p h 
respectively, as well as subscripts + or — for their symmetric and anti-symmetric 
combinations, i. e., ip± = ipi ± ip v = ip(pi) ± ip{p v ) for any quantity ip. 

By differentiation of the free energy we get analogous relations for the pressure P 
and chemical potential p, the two quantities that directly enter the determination 
of the coexistence curve: 

P(p + Ap, a u + Ao"h) = P{p) + A^h P{p) + Ap P\p) + ..., 



19 



H(p + Ap, tr H + A<r H ) = A*(p) + Aa H fi(p) + Aa H p'(p) + . . . . 

For given an, the coexisting densities are obtained from the criterion of equal 
pressure and chemical potential, P(p v ,an) = P(pi,an) and p(p„, an) = p(pi,an)- 
Specializing to the coexisting densities p v and pi obtained for this an, the equilib- 
rium conditions reduce to P v = Pi and p v = pi, or 

p_ = /1 _= o. (Al) 

Non-zero Aan, on the other hand, modifies both pressure and chemical potential 
and is therefore generally accompanied by changes Ap„ and Ap^ in the coexist- 
ing densities. Taking Eq. (|A1J) into account, switching to symmetric and anti- 
symmetric combinations of quantities, and simplifying, the equilibrium conditions 
are easily obtained as 

PL Ap+ + P' + Ap_ + 2 P_ Act h = 0, 

fj!_ Ap + + fjf + Ap_ + 2 (i- Aa H = 

to lowest order. Solving for Ap± yields the sensitivity of the symmetric and asym- 
metric shifts in the coexisting densities to a change Aan of the effective diameter 
as 

a p± 2 pLp--p' ± P- 

Aan PLp'±~PLp' t 
P' ± P--p' ± P- 
pi ,,i — pi ,,i ' 

where we have taken advantage of some cancellations in the denominator to arrive 
at the second expression. 

From the product of two density derivatives in the denominator as opposed to 
single factors of density derivatives in the numerator it is already apparent that the 
sensitivity of the binodal to variation of an must be proportional to the isothermal 
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compressibilities k v and Ki at the coexisting densities. To make this more explicit 
and simplify the expressions further, we introduce the free energy per unit volume 
/ and use basic thermodynamic relations to replace the pressure and chemical 
potential, as well as their derivatives, by appropriate expressions in terms of / and 
its derivatives while converting reciprocals of /" to the isothermal compressibilities. 
After some slightly tedious simplification we get explicit expressions for the shifts 
in the coexisting densities induced by a change A<th of the effective diameter, viz., 



The change of any one of the coexisting densities with o"h is thus proportional to 
the product of the compressibility at that density and a factor that combines the 
sensitivities of the free energy at both branches of the binodal with its density 
derivative. 

In an exact theory, of course, all the <th derivatives vanish exactly and the 
coexisting densities are strictly independent of <7h even at the critical point, (3 = j3 c , 
P- = 0, where the compressibility diverges. In an approximate theory, on the other 
hand, /' — /_/ Ap_ is non-zero for (3 > j3 c , and while its limit for (3 — > (3 C obviously 
vanishes, that of its product with the diverging compressibility may be zero, finite, 
or even infinite; for BAS SCOZA, Figs. El and El indicate a finite value. As (3 further 
increases from j3 c , the two branches of the binodal separate rapidly, and the one- 
sided difference quotient /_/Ap_ over the coexistence region no longer cancels 
the derivative /' at either side of the binodal. The expression in parentheses is 
therefore no longer expected to be close to zero. In the case of BAS SCOZA, 
the variability of / with <th strongly increases with density, as is expected for 
thermodynamic perturbation theories in general. This carries over to /' — /_/Ap_ 
and, in combination with the prefactor p 2 , so explains why the sensitivity of p\ is 




Ap; Ap„ 

Aa H Aa H ' 
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Tables 



Tab.|TJ Sensitivity of free energy, pressure, and chemical potential at p* — 1 and 
p* = 0.7 to the choice of <th- For any quantity X, AX is the difference between 
the results obtained with obfc — 2.3 and obfc — 2.5. The columns AX/X and 
AX/ (dX/ dp) relate the response of both the value and the slope of X as a function 
of p and so characterize the sensitivity of the shape of X(p) to variation of obfc- 
- In general, the free energy is seen to be far more forgiving than P and p with 
respect to the precise choice of o"h, especially at low temperature. 
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Figures 



Fig. Dependence of the BAS SCOZA excess free energy v4 cx on the effective 
hard core diameter <th for the four states considered in Ref. |l6| . (This corresponds 
to Fig. 2 of Ref. 0.) 

Fig. 121 The pressure P as a function of p along four isotherms as computed 
within BAS SCOZA with the BFC diameters with obfc — 2.3, 2.5 and e, as well 
as the simulation data of Ref. 2^1- (This corresponds to Fig. 5 of Ref. The 
inset gives a detailed view of the isotherm at T* = 1.4, along with the curves 
obtained in an MSA based variant of the theory (K = (3) for czbfc — 2, 2.3, and e. 

Fig. The excess free energy A cx as a function of p along four isotherms as 
computed with the BFC diameter with gibfc — 2.3, 2.5, and e. 

Fig. EJ Dependence of the chemical potential u on the effective hard core di- 
ameter cth for the four states considered in Ref. [16j ]. 

Fig. EJ The curve of phase coexistence as obtained for czbfc = 2.2, 2.3, 2.5, and 
e as well as with the BH diameter, along with simulation results of Ref. y. ah 
data refer to a truncated but not shifted, rather than the full, LJ potential. 

Fig- [HI The curve of phase coexistence as obtained for gibfc — 2.2, 2.3, 2.5, and 
e as well as with the BH diameter, along with simulation results of Ref. [26] . All 
data refer to the full LJ potential. 

Fig. [3 SCOZA (Af A ) and BAS (AA S ) contributions to the BAS SCOZA excess 
free energy at T* = 1.35, p* = 0.65 for varying <th- Clearly, the SA free energy 
A ex = Aff A + A As is far less sensitive to cth than either A^ A or A As taken by 
itself. On the other hand, A c * is not quite constant as must be the case for an 
exact theory. 
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Fig. |H| Diameter dependence of the internal energy at four different states. 



Fig. EJ SCOZA (U^ A ) and BAS (AU S ) contributions to the BAS SCOZA in- 
ternal energy at T* = 1.35, p* = 0.65 and T* = 1.5, p* = 0.3 for varying a#. 



Fig. E3 BAS SCOZA and simulation |25j results for the configurational in- 
ternal energy U cx along four different isot 
corresponds to one panel of Fig. 4 of Ref. 



lerms for three choices of obfc- (This 

m 
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p* T* 


1 AA ex 1 AA ox 


1 AP 1 AP 


A^ Ap 


4° dA c */dp* 


1 P 1 dP/dp* 


p dp/ dp* 


0.7 4.0 
2.5 
1.4 
0.75 


2.6 • 10~ 3 7.8 • 10" 4 
5.0 • 10~ 3 5.2 • lO" 4 
8.5 • 10~ 4 2.5 • 10~ 3 
6.7- 10~ 4 7.6- 10" 4 


8.0 • 10" 4 2.0 • 10" 4 

1.1 ■ 10~ 3 2.3 • 10" 4 

1.2 • 10~ 2 1.4 • 10- 3 
3.5 • lO- 2 2.6 • 10" 2 


9.5 • 10~ 3 1.7 • 10" 3 
1.0 • lO" 2 1.0 • 10~ 3 

2.6 • 10~ 3 3.1 • 10~ 4 
6.6 • lO" 3 2.3 • lO" 2 


1.0 4.0 
2.5 
1.4 
0.75 


4.0- 10" 4 1.4- 10~ 4 
3.2 • 10~ 3 7.1 • 10" 4 
6.2 • 10~ 2 2.4 • 10- 3 
8.0 • 10~ 3 6.2 • 10" 3 


4.8 • 10~ 3 1.2 • 10- 3 

1.0 • ID -2 2.2 • lO -3 

2.1 • lO" 2 3.5 • 10~ 3 
4.4 • lO" 2 4.6 • lO" 3 

TABLE I: 


2.1 • 10~ 3 5.9 • 10" 4 
8.3 • 10~ 3 1.9 • 10~ 3 
2.5-10-2 3.6- 10~ 3 
2.4- lO" 1 5.0- 10~ 3 
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